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Abstract. We give an overview of the progress that has been made in recent years 
in understanding the dynamic multiscaling of homogeneous, isotropic turbulence and 
related problems. We emphasise the similarity of this problem with the dynamic scaling 
of time-dependent correlation functions in the vicinity of a critical point in, e.g., a spin 
system. The universality of dynamic-multiscaling exponents in fluid turbulence is 
explored by detailed simulations of the GOY shell model for fluid turbulence. 



PACS numbers: 47.27.Gs, 47.53.Eq 
1. Introduction 

The statistical properties of fully developed, homogeneous, isotropic turbulence are often 
characterised by velocity structure functions, which are averages of the differences of 
fluid velocities at two points separated by a distance r (a precise definition is given 
below). If r lies in the inertial range of scales that lie between the large length scale 
L, at which energy is pumped into a turbulent fluid, and the small dissipation scale 
rjd, at which viscous dissipation becomes significant, these structure functions scale as 
a power of r, in a manner that is reminiscent of the algebraic behaviour of correlation 
functions at a critical point in, say, a spin system. This similarity between the statistical 
properties of homogeneous, isotropic turbulence and a system at a critical point is well 
known; and its elucidation for equal-time structure functions has been the subject of 
many papers: It turns out that the simple scaling we are accustomed to at most critical 
points must be generalised to multiscaling in turbulence; i.e., an infinity of exponents is 
required to characterise the inertial-range behaviours of structure functions of different 
orders []]. 

| Author to whom all correspondence should be addressed 
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The scaling behaviour of correlation functions at a critical point [2J is associated 
with the divergence of a correlation length £ at the critical point; e.g., in a spin system 
£ ~ t v if the field H = 0, where the reduced temperature t = (T — T c )/T c , T is the 
temperature, T c the critical temperature, and z/ is an equal-time critical exponent that is 
universal (for systems in a given universality class). In the vicinity of a critical point the 
relaxation time r, which can be determined from time- dependent correlation functions, 
scales as follows: 

r~£*- (1) 

This is known as the dynamic-scaling Ansatz via which we define z, the dynamic-scaling 
exponent. Over the past few years there has been considerable progress [3j HI [5J El El El [9] 
in developing the analogue of such a dynamic-scaling Ansatz for izme- dependent structure 
functions in homogeneous, isotropic turbulence. We give an overview of this work on 
dynamic mutiscaling in turbulence. We show, in particular, that (a) an infinity of 
dynamic-multiscaling exponents is required here, (b) these exponents depend on the 
precise way in which relaxation times are extracted from time-dependent structure 
functions, and (c) that dynamic-multiscaling exponents are related by linear bridge 
relations to equal-time multiscaling exponents. 

The remaining part of this paper is organised as follows: Section 2 contains an 
introduction to the multiscaling of equal-time structure functions in fluid turbulence. 
Section 3 is devoted to the dynamic multiscaling of time-dependent structure functions. 
In Sect. 4 we introduce the GOY shell model pQ [101 [11] for fluid turbulence and 
give the details of our numerical studies of this model. Section 5 gives representative 
results from our simulations for time-dependent structure functions for the GOY shell 
model with conventional viscosity; we also present new results with hyperviscosity and 
show explicitly that dynamic multiscaling exponents are independent of the type of 
viscosity we use. We end with concluding remarks about the possibility of experimental 
verifications of our predictions and generalizations to other types of turbulence such as 
passive-scalar turbulence [12] . 

2. Equal-time Multiscaling 

Fluid flows are described by the Navier-Stokes equation for the velocity field u(x, t) at 
point x and time t: 

9 ( u + u.Vu = -VP + i/ V 2 u + f; (2) 

we consider low-Mach number flows that are nearly incompressible and so equation ([21) 
must be augmented by the incompressibility constraint 

V.u = 0, (3) 

which can be used to eliminate the pressure P in equation (j2]); we choose the uniform 
density p = 1; z/q is the kinematic viscosity; and for decaying turbulence the external 



Dynamic Multiscaling in Turbulence 



3 



force f is zero. Turbulence occurs when the Reynolds number Re = (£v)/u is large; 
here I and v are characteristic length and velocity scales of the flow. 

The chaotic nature of turbulence has led to a natural statistical description of 
the velocity field in terms of velocity structure functions. For instance, we can 
define the order-p, equal-time, structure function for longitudinal velocity increments 
5u\\ (x, r, t) = [u(x + r, t) — u(x, t)} .r/r as follows: 



the power-law dependence on r, which defines the exponent ( p , holds for L <C r <C T]d] 
and the angular brackets indicate an average over either the statistical steady state, for 
forced turbulence, or statistically independent initial conditions, for decaying turbulence. 

In contrast with the scaling behaviour of correlation functions at conventional 
critical points in equilibrium statistical mechanics, in turbulence the structure functions 
S p (r) do not exhibit simple scaling forms: Experimental and numerical evidence suggests 
that S p (r) show multiscaling, with ( p a nonlinear, convex, monotone increasing function 
of p pp. The 1941 theory of Kolmogorov (K41) [131 Ej yields simple scaling with 
Cjf 41 = p/3; the measured values of ( p deviate significantly from Cjf" 41 f° r V > 3; and for 
p = 3 we have the exact von Karman-Howarth [1] result (3 = 1. 

Even though the K41 phenomenology fails to capture the multiscaling of the equal- 
time structure function S p (r), it provides us with important conceptual underpinnings 
for studies of homogeneous, isotropic turbulence. In particular, the exponents C^ 41 are 
universal in the sense that they do not depend on the details of the dissipation, i.e., 
the viscosity. (Of course measurements should be made far away from boundaries so 
that the flow satisfies the conditions of homogeneity and isotropy.) This suggestion 
of power-law scaling with universal scaling exponents was made for turbulence a few 
decades before it was appreciated fully in the context of critical phenomena. 

The universality of the exponents ( p holds even if we go beyond K41 phenomenology. 
Let us first examine the dependence of these exponents on the dissipation mechanism. 
In numerical simulations it is possible to introduce a hyperviscosity v a by replacing the 
viscous term u V 2 u in equation ([2]) by z/ a V 2+a u; here a > determines the degree of 
hyperviscosity; and a = yields normal viscous dissipation. Some early shell-model 
studies [151 EEl ELY] suggested that the exponents ( p depend on a. But subsequent direct 
numerical simulations (DNS) [H] of equations (EH) and numerical studies of a shell 
model [19] for turbulence have argued against this. Our results are consistent with 
these later studies. 

The exponents ( p do not seem to depend on whether they are measured in 
statistically steady, forced turbulence or in decaying turbulence. Numerical evidence 
for this universality has been provided by the shell-model studies of Ref. [20] . 

In the remaining Sections of this paper we will show how to generalise the dynamic- 
scaling Ansatz (OQ) to account for multiscaling in turbulence. Our discussion will be based 
on earlier studies [3], [6] and the work carried out in our group [HE]. We will then examine 
the universality of the dynamic-multiscaling exponents in fluid turbulence, i.e., we will 
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explore their dependence on (a) the hyperviscosity parameter a and (b) the type of 
turbulence (statistically steady as opposed to decaying). 

3. Dynamic Multiscaling 

The dynamic scaling of time- dependent correlation functions at a critical point in an 
equilibrium system was systematised soon after the scaling of equal-time correlation 
functions. The analogous development for homogeneous, isotropic turbulence has been 
carried out recently [31 HI El El El- We summarise the essential points here before 
presenting our new results. A naive extension of K41 phenomenology to time-dependent 
structure functions yields a dynamic exponent z KA1 = 2/3 for all orders p. Two 
improvements are required to go beyond this K41 result: (a) We must account for the 
multiscaling of velocity structure functions, (b) We must distinguish between structure 
functions of Eulerian and Lagrangian velocities; the former yield trivial dynamic scaling 
with z £ = 1, for all p, since the mean flow (or the flow of the largest eddy) directly 
advects small eddies and so temporal and spatial separations are related linearly by 
the mean flow velocity; nontrivial dynamic multiscaling can be anticipated, therefore, 
only for Lagrangian [21] or quasi-Lagrangian [6] velocity structure functions. The latter 
are defined in terms of the quasi-Lagrangian velocity u that is related to its Eulerian 
counterpart u as follows: 

u(x,t) = u[x + R(t;r o ,0),*], (5) 

with R(i; ro, 0) the position at time t of a Largrangian particle that was at ro at time 
t = 0. Equal-time, quasi-Lagrangian velocity structure functions are the same as their 
Eulerian counterparts [22J. 

The order-p, time-dependent, structure function, for longitudinal, quasi-Lagrangian 
velocity increments is [31 H] 

F p (r, {ti, . . . ,t p }) = ([<ytt||(x,ti,r) . . . 5u\\(x, t p , r)]) . (6) 

Since we are interested in scaling behaviours, we restrict r to the inertial range, consider, 
for simplicity, t\ — t and ti = . . . = t p = 0, and denote the structure function (jHJ) by 
J- p (r, t). Given !F p {r : t), there are different ways of extracting time scales. For example, 
we can use the following order-p, degree-M integral- and derivative-time scales that are 
defined, respectively, as [1] 



and 



r p D M {r) 



1 



Sp(r) 



F p (r,t)t {M - l) dt (7) 



1 d M T p {r,t) 



S p (r) dt M 



t=oJ 

If the integral in (JTj) and the derivative in (jHJ) exist we can generalise the dynamic- 
scaling Ansatz §B) at a critical point to the following dynamic-multiscaling Ansatze for 
homogeneous, isotropic turbulence 0, H] : 

?mW~^ m ; (9) 
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r p D M {r) ~ r**. (10) 

These equations define, respectively, the integral- and the derivative-time multiscaling 
exponents z* M and z® M , which satisfy the bridge relations 

4,M = 1 + [Cp-M ~ Cp]/M, (11) 

obtained first in Ref. [3], and 

z?,m = 1 + [Cp-Cp+m]/M, (12) 

obtained in Ref. [3|. These bridge relations follow from a generalisation [HI H] of the 
multifractal formalism [1] for turbulence. Here we sketch the arguments that lead to 
equation (TTTI) for statistically steady turbulence and refer the reader to Refs. [11 [23] 
for details: We begin by assuming the following multifractal form for the order-p time- 
dependent structure functions: 

oc f d^h)C-f +ph - v ^GV' h {—). (13) 

U L Jl L T P,h 

As in the equal-time multifractal formalism [TJ, the scaling exponents h e X = 
[hmin, h max }. Corresponding to each exponent h there is a set C M 3 , of fractal 
dimension D u ( h > and with a measure dfi(h), such that 5u r (x.)/uL ~ (r/L) h if x 6 S^, 
with 5u r (x) = |u(x + r) — u(x)|. We assume furthermore that the scaling function 
gp,h^_t_^ j g sucn ^ na ^ Q p ' h (0) = 1 and that the characteristic decay time T Pj h ~ 

r/6u r (x) ~ r 1_/l . If we substitute for jF p (r, t) in equation flU]), do the time integral first 
by a saddle-point method, we obtain the bridge relation ffTTj) . Similar calculations lead 
to (}I2"]) . A complete discussion of such bridge relations, as well as similar relations for 
passive-scalar turbulence, can be found in Refs. [U [23] • The last of these references 
also shows that the bridge relations are the same for decaying and statistically steady 
turbulence. 

Given current computational resources, it has not been possible to verify the 
bridge relations ( fTHfT2l) by computing quasi-Lagrangian velocities in a direct numerical 
simulation of the Navier-Stokes equation (J2H3D f° r an incompressible fluid. Thus we 
must turn to numerical studies of simple shell models of turbulence that we describe in 
the next Section. 

4. The GOY Shell Model 

We have carried out extensive numerical simulations to obtain equal-time and dynamic 
multiscaling exponents for the GOY shell model for fluid turbulence [U [101 EH] • This 
shell model is defined on a logarithmically discretized Fourier space labelled by scalar 
wave vectors k n that are associated with the shells n. The dynamical variables are 
the complex scalar shell velocities u n (k n ), henceforth denoted by u n . The evolution 
equations for the GOY model are 



4 + e 



dt \ k. 



Ur. 



a n U n+ iU n+2 + + C n M„_iM r 



(14) 
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where k n = k 2 n , k = 1/16, complex conjugation is denoted by *, and the coefficients 
a n = k n , b n = —5k n -i, c n = —(1 — 5)k n _ 2 are chosen to conserve the shell- model 
analogues of energy and helicity in the inviscid, unforced limit, and 1 < n < N . We 
use the standard value 5 = 1/2. The second term allows for hyperviscous dissipation 
with degree a (the case a = corresponds to conventional viscous dissipation); kd is a 
large wavenumber whose inverse is comparable to the dissipation length scale; we choose 
kd = 2 18 and a = or a = 2. Since we concentrate on decaying turbulence here, the 
forcing term, required to drive the system into a statistically steady state, is absent. 
The logarithmic discretization of Fourier space allows us to reach very high Reynolds 
numbers in numerical simulations of the GOY model even with N = 22. 

We use a slaved, second-order, Adams-Bashforth scheme [2JJ [25] to integrate the 
GOY-model equations with N = 22 shells, and, in some representative cases, N = 35 
(this requires a fifth-order scheme [26J ), with the boundary conditions u n = for n < 
and n > N. In our simulations we use 5t = 10 -4 as the integration time step, and the 
viscosity v = 10~ 7 for a = and v = 10~ 6 for a = 2. 

The GOY-model equations allow for direct interactions only between nearest- and 
next-nearest-neighbour shells. By contrast, in the Fourier transform of the Navier- 
Stokes equation every Fourier mode of the velocity is coupled to every other Fourier 
mode directly. This direct sweeping, as it is called, leads to the trivial dynamic scaling 
of Eulerian velocity structure functions (Sect. 3). Since the GOY model does not have 
direct sweeping in this sense, it is sometimes thought of as a highly simplified quasi- 
Lagrangian version of the Navier-Stokes equation. Thus we might expect nontrivial 
dynamic multiscaling for GOY-model structure functions. The equal-time structure 
functions for this model are 

S p (k n ) = ({u n (t)u* n (t)f 2 )~k^, (15) 

where the power-law dependence is obtained only if fc" 1 lies in the inertial range. 
Three cycles [27] in the static solutions of the GOY model lead to rough, period-three 
oscillations in S p (k n ); thus we use the modified structure function [27] 

S p = (\%[u n+2 u n+1 u n - (l/4)n n _iu n u„ +1 ]| p/3 ) 

~ K (l % (16) 

in which these oscillations are effectively filtered out. Therefore we use equations ffT6l) 
and ( ITTj) to extract ( p . Time scales are obtained from the order-p, time-dependent 
structure functions for the GOY model, namely, 

F p {k n ,t Q ,t) = (K(t R(t + t)F /2 ). (17) 

In Refs. [U [23] we have used such time- dependent structure functions to verify 
the bridge relations fllltil2p for statistically steady turbulence. Here we give a short 
overview of our results for decaying turbulence with conventional viscosity (a = 0) and 
hyperviscosity (a = 2). In our studies of decaying turbulence we have used two types 
of initial conditions; in both of these all the energy is initially concentrated in the first 
few Fourier modes, i.e., at large length scales, as in typical wind-tunnel experiments of 
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homogeneous, isotropic turbulence pp. In the first initial condition u n = kn 2 e ien , for 
n = 1, 2, and u n = kl/ 2 e~ k " 2 e l0n for 3 < n < N, where 8 n is a random phase angle 
distributed uniformly between and 2tt; we use this for the case a = 0. The second 
initial condition we use is similar, namely, u n = e~ k ™e l9 " for 1 < n < N . We have used 
both these initial conditions for the case a = 2; our results do not depend significantly 
on the initial condition that is used; the results we report here for a = 2 are for the 
second type of initial condition. 

Time is measured in terms of the initial large eddy-turnover time t^. For the GOY 
shell model = l/(u rms ki), where the root-mean-square velocity u rms = [(%2 n l^nl 2 )] 1 ^ 2 
is defined for the initial velocity field. The shell-model energy spectrum is [23] 
E{k n ) = (|w n | 2 ) /k n = S2(k n )/k n ; since S^t^n) shows the period-three oscillations 
mentioned above, a smooth energy spectrum is obtained by using T, 2 {k n )/k n as we show 
below. The mean kinetic energy dissipation rate and the integral scale for the GOY 

model are, respectively, e = v{k n /k d ) a kl\u 2 n \) and L int = 7^7^77^71 • As the 

turbulence decays, Li nt increases; once it becomes comparable to the size of the system 
the total energy decays [H] as t~ 2 . Our results for dynamic- multiscaling exponents 
are obtained for times that are much shorter than the time over which L» n t becomes 
comparable to the system size (~ k^ 1 for the GOY model). 

5. Results and Conclusions 

We begin by looking at the mean energy dissipation rate e as a function of time. A 
representative plot, averaged over 2000 initial conditions, is shown for a = 2 in Figure 
[TJ Though this plot is noisy, it shows a clear peak. This peak, at t/t L ~ 1.2 (t L ~ 5) 
in Figure [D, signals the completion of the cascade that transfers energy from the scale 
at which it is injected to the small scales where viscous dissipation becomes significant. 
Representative energy spectra at cascade completion are shown in Figure [2] after this 
point in time the energy spectrum decays very slowly without an appreciable change in 
the slope of the scaling regime, which appears as a nearly straight-line segment in the 
log-log plots of Figure EJ Straight-line segments also show up in plots of the structure 
function S p at (or after) cascade completion as shown in the representative plots, for 
a = 2, of Figure[3](a); each curve in this Figure has been averaged over 5000 independent 
initial conditions. From the slopes of the straight-line segments in these plots (see Table 
1) we obtain the equal-time multiscaling exponents ( p whose dependence on p is shown 
in Figure [3](b). The values we quote for these (and other) exponents are the means of 
the slopes of 50 different plots like Figure [3](a), which are obtained from 50 statistically 
independent runs; the corresponding standard deviations yield the error bars shown in 
Table 1. 

Our results for the equal-time exponents ( p for a = and a = 2 are presented 
in Table 1. By comparing Columns 2 and 3 in this Table we see that the exponents 
for both values of a agree with each other and with the exponents reported earlier for 
statistically steady turbulence [3j. Thus we reconfirm the universality of equal-time 
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Figure 1. The mean energy dissipation rate e versus time for a — 2 (for clarity 
we show data for 0.8 < t/t^ < 2.5); these data have been averaged over 2000 initial 
conditions. The main peak at tjth — 1-2 is a signature of cascade completion. For 
a = e displays a similar peak. 




Figure 2. (a) Log-log plot of the kinetic energy spectrum E(k n ) = S2(k n )/k n versus 
k n with the period-three oscillations (see text). The black dashed line corresponds 
to the K41 scaling prediction E{k n ) ~ k~ 5 ^ 3 . (b) Log-log plot of Y,2{kn)/k n versus 
fc„; note that the period-three oscillations are suppressed here; the black dashed line 
indicates the K41 scaling prediction. 



exponents: they neither depend on the precise dissipation mechanism nor on whether 
we consider statistically steady or decaying turbulence. 

In Ref. [1] the dynamic-multiscaling exponents z 1 M and z® M were obtained, by 
using integral- and derivative-time scales, for statistically steady turbulence in the 
GOY model with a = 0. We give below an overview of our recent results for dynamic 
multiscaling but for decaying turbulence; for a detailed discussion of these results we 
refer the reader to Ref. |23j . 

Time- dependent structure functions in decaying turbulence must, of course, depend 
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Figure 3. (a) Log-log plot of the modified structure function £ p , for p = 1 to 4 (top 
to bottom), versus k n for a = 2 and number of shells N = 35 (we show data only 
for the first 22 shells), (b) Plot of the equal-time multiscaling exponents ( p (obtained 
from (a) and listed in Table 1) versus the order p. The open circles (o), connected by 
a curve to guide the eye, indicate data from our numerical simulations; the thick black 
line is the K41 prediction Cf 41 = p/3. 



Table 1. Our results for the equal-time multiscaling exponents Cp a t cascade 
completion for a = (Column 2) and a = 2 (Column 3). We indicate, in parentheses, 
the ranges of shell numbers n over which we fit our data for equal-time structure 
functions to obtain these exponents. Note that the exponents in Columns 2 and 3 
agree very well with each other. 



order-p 


Cv (4-14) 


(p (4-16) 


1 


0.380 ± 0.001 


0.37 ± 0.01 


2 


0.709 ± 0.003 


0.699 ± 0.008 


3 


1.000 ± 0.005 


1.003 ± 0.008 


4 


1.266 ± 0.008 


1.29 ± 0.02 


5 


1.51 ± 0.01 


1.55 ± 0.03 


6 


1.74 ± 0.02 


1.79 ± 0.05 



on the origin of time to from which we start making measurements. It turns out that this 
dependence can be eliminated by normalising these structure functions by their values 
at the origin of time; we have shown this analytically for the Kraichnan model [281 EE] 
for passive-scalar turbulence and numerically in several other models [23]. We present 
some representative numerical results for the GOY model below for which we use the 
normalised time-dependent structure functions 

F p (k n ,t ,t) 



Qp(k n , t) 



F p (k n ,t ,0)' 



It turns out that Q p does not depend on t as shown explicitly in Figure @] for p = 4 



Dynamic Multiscaling in Turbulence 



10 



1 

0.9 

? - 8 
0*0.7 

0.6 

O.5I ' 1 

0.1 0.2 0.3 

Figure 4. Plots of Qi(kn,t/tL) for n = 5 obtained from 6 different values of to] 
successive values of to are separated by 0.5£l- Since the symbols used for different 
values of to are indistinguishable on the scale of this figure, we conclude that Q4 does 
not depend on to- 

and n = 5. This is why we have not displayed to as an argument of Q p . 

Given that Q p does not depend on to, we can use expressions such as (l9l HD]L with T p 
replaced by Q p , to obtain dynamic-multiscaling exponents for decaying turbulence. We 
consider here the order-1 integral- and order-2 derivative-time scales and the associated 
exponents z pl and z p2 , respectively. 

To calculate the integral-time scale T pl , we evaluate the integral in equation (JTj) 
with Q p instead of T p and the upper limit replaced by t^, the time at which Q p (k n ,t) 
reaches a value /i G [0, 1]. We should, of course, use /J — 0, i.e., = 00, but it is hard 
to obtain well-averaged data for large t, so we use \i = 0.6. We have checked that our 
results are not affected significantly if we restrict ourselves to the range 0.4 < // < 0.7. In 
Figure 0(a) we show representative plots of Q 4 versus time t/t L for shells n = 3,4, 7, and 
9 and a = 2; from this we obtain T/ x as described above. The slopes of log-log plots of 
Tp i versus k n now yield z pl as shown in Figure Mb) for p — 4. The derivative-time scale 
T^ 2 is calculated by using a sixth-order, centred difference scheme. The derivative-time 
exponent z p2 is then obtained from slopes of log-log plots of T p2 versus k n . 

Tables 2 and 3 summarise our results for the dynamic-multiscaling exponents. 
All these exponents have been calculated for k n in the inertial range 4 < n < 14; 
the exponents z p2 are shown only for a = 0. In Tables 2 and 3 we also give the 
exponents that we get by using the bridge relations ( fTTl) and ( fl~2l) along with the equal- 
time exponents given in Table 1; there is good agreement between the exponents from 
our numerical simulations and those obtained via the bridge relations. Moreover, the 
exponents z pl and z® 2 , for a = 0, and z pl , for a = 2, which we obtain for decaying 
turbulence, are equal (within error bars) to their counterparts (see Table 2 of Ref. [1]) 
for statistically steady turbulence. Thus the dynamic-multiscaling exponents for the 
GOY model are universal. Plots of these exponents versus p are shown in Figures HD^a) 
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Figure 5. (a) Illustrative plots of the normalised, fourth-order, time-dependent 
structure function Qi(k n ,t) versus time tjt^ for shells n = 3,4,7 and 9 (from top 
to bottom), (b) A log-log plot of the integral-time scale T\ x versus k n \ the integral- 
time exponent z\ 1 — 0.728 ± 0.006 follows from the slope of the line (in the range 
4 < n < 12). 



Table 2. Column 2: The dynamic multiscaling exponents z£ l from the bridge relation 
(lll[) and our numerical results for £ p (Table 1). Column 3: The same exponents as in 
Column 2 but now obtained from our numerical simulations of decaying turbulence in 
the GOY shell model (for the case a — 2). Note the agreement between the exponents 
in Columns 2 and 3. 



order-p 


^[Eq.tflT])] 


z 1 


1 


0.630 ± 0.001 


0.62 ± 0.01 


2 


0.675 ± 0.009 


0.683 ± 0.003 


3 


0.69 ± 0.01 


0.712 ± 0.004 


4 


0.71 ± 0.02 


0.728 ± 0.006 


5 


0.74 ± 0.05 


0.750 ± 0.009 


6 


0.76 ± 0.08 


0.76 ± 0.01 



and [6](b). 

The universality of dynamic exponents also goes through for models of passive- 
scalar turbulence as we show in Refs. [5j [23]. It turns out that the Kraichnan 
model [2HI [29] shows simple scaling. However, a shell-model version of a passive 
scalar advected by GOY-model velocities exhibits dynamic multiscaling with dynamic 
exponents that depend on the degree M, defined in equations ([7}[8]), but not on the 
order p. 

We hope our detailed numerical simulations of dynamic multiscaling in the GOY 
shell model for fluid turbulence will lead to experiments designed to measure time- 
dependent Lagrangian or quasi-Lagrangian structure functions in fluid turbulence. 



Dynamic Multiscaling in Turbulence 



12 



Table 3. Dynamic multiscaling exponents from our simulations of decaying turbulence 
in the GOY model with a = 0: The bridge relations (|imi2p are used along with the 
equal-time exponents Cp for a = (Table 1) to calculate 1 (Column 2) and z® 2 
(Column 4). The corresponding dynamic- multiscaling exponents that we obtain from 
our numerical study of the normalised, time-dependent structure functions Q p are 
given in Columns 3 and 5. 



order-p ^ jEq.(fl~Tj)] 


z 1 
p,i 


z° 2 {Eq.m\ 


z D 


1 0.620 ± 0.001 


0.60 ± 0.02 


0.690 ± 0.006 


0.687 ± 0.003 


2 0.671 ± 0.004 


0.67 ± 0.03 


0.72 ± 0.01 


0.719 ± 0.005 


3 0.709 ± 0.008 


0.707 ± 0.006 


0.74 ± 0.02 


0.743 ± 0.007 


4 0.73 ± 0.01 


0.736 ± 0.008 


0.76 ± 0.03 


0.75 ± 0.01 


5 0.75 ± 0.02 


0.752 ± 0.009 


0.77 ± 0.03 
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Figure 6. Plots of the dynamic multiscaling exponents (a) 1 and (b) z® 2 versus 
p; the error bars are obtained as described in the text. These plots compare data for 
statistically steady turbulence in the GOY model (from Table 2 in Ref. [4]) and our 
data for decaying turbulence in the GOY model with normal viscosity (a — 0) and 
hyperviscosity (a — 2); and they illustrate the universality of dynamic-multiscaling 
exponents that we discus in the text. 



Advances in experimental techniques have made it possible to get good data for, say, 
the acceleration of Lagrangian particles [30] . We believe such techniques can be refined 
to measure the sorts of time-dependent structure functions we have discussed here. 
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